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The parity doublet model, containing the SU(2) multiplets including the baryons iden- 
tified as the chiral partners of the nucleons is applied to neutron stars. The maximum 
mass for the star is calculated for different stages of the cooling taking into account 
finite temperature/entropy effect, trapped neutrinos and fixed baryon number. Rotation 
effects are also included. 
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1. Introduction 

It has been shown that effective chiral hadronic models can satisfactorily describe 
nuclear matter, properties of finite nuclei as well as the structure of neutron starJl^. 
While most of the approaches are based on variations of the linear sigma model or 
its non-linear realization, in the approach discussed here the scalar sigma meson 
serves to split the nucleon and its chiral partner (particle with opposite parity), 
while in the chirally restored phase both baryonic states become degenerate but 
not massless. The presence of the chiral partner allows to have a bare mass term 
mo in the Lagrangian density in such a way that it does not break chirality because 
the physical fi elds in this case are a mixture of the fields of the particles and their 
chiral partner 
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2. The parity Doublet Model 

It is assumed that the star is in chemical equilibrium and the baryons interact 
through the mesons cr, u> and p (included in order to reproduce the high asymmetry 
between neutrons and protons). Electrons are included to insure charge neutrality. 
The Lagrangian density in mean field approximation becomes 

LmFT = Lkin + Lb seal + Lbvcc + -^tiec + Lg^al + -^SS, (1) 
Lsscal + Lbvcc = - ^ i^ddiujlOi^^ + giplQTsP^ + M*]ll:i, (2) 



= -hmlu:^ + my) - g,[u^ + Gp^c^ + p% (3) 



2 

Lscal = ^M^fT^ - (4) 

LsB = mlf^a, 

where besides the kinetic term for the fcrmions there are terms of interaction be- 
tween the baryons and the scalar and vector mesons, self- interaction terms for the 
vector and scalar mesons and an explicit symmetry breaking term included in order 
to reproduce the masses of the pseudo-scalar mesons. 

The effective masses of the nucleons and their chiral partners reproduce their 
measured values at low densities, when the scalar condensate has its vacuum value, 
and they go to a specific value mo = 790 MeV (chosen in order to have a physical 
compressibility at saturation"^) at high densities, when the scalar condensate go to 
zero 



Ml = 



M,^. (5) 

2 (To 



A possible, but not definite, candidate for the nucleon chiral partner is the 
A^_(1535), but besides that, the case with A^_(1200) is included to study the effect 
of the mass of the N_. This variation has drastic consequences on the results. For 
this model, four different cases first studied in SI are applied to proto-neutron and 
neutron stars: 

• PI: Mn_ = 1200MeV and 54 = 

• P2: Mn_ = 1200MeV and 54 = 200 

• P3: Mn_ = 1500MeV and 54 = 

• P4: Mn_ = ISOOMeV and 54 = 200 



As will be shown, the results are also very sensitive to the value of the coupling 
constant of the self-interaction of the vector mesons (74. 
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Fig. 1. Scalar condensate as a function of chemical potential for different configurations 

3. Neutron Stars 

For the parity model, when the density increases, the protons, the neutron and 
the proton chiral partners appear, respectively. Since the matter is not symmetric, 
in general the two isospin states of the chiral partner appear at different densities, 
which induces the chiral restoration to happen earlier compared with the symmetric 
nuclear matter case. As can be seen from FigUthe scalar condensate used in this 
case as the order parameter for the chiral restoration shows that this transition 
happens very smoothly, and can be identified as a crossover for all cases. This effect 
is caused by the introduction of beta equilibrium and charge neutrality. It can also 
be seen in this plot that for P3 the chiral restoration caused by the appearance of 
the chiral partners happens at very high densities, beyond the ones possible inside 
a neutron star. 

The maximum mass of the star is calculated solving the TOV equations and it 
shows that higher maximum values are reached for 54 equal to zero (Figl2|). That 
means that when (74 increases, the value of the vector meson co related to it decreases 
its value and consequently its repulsive effect, in such a way that the star can hold 
smaller quantity of mass against collapse. Besides the EOS for the dense part of 
the star, the EOS for an outer crust, an inner crust and an atmosphere have been 
addecP. 

4. Cooling 

To study the cooling of such complex system two different features are taken into 
account separately: finite entropy per baryon and lepton number. The reason for 
considering finite entropy instead of temperature comes from the fact that in this 
case the temperature is higher in the center of the star and colder at the edge, which 
is more realistic than assuming the whole star at equal temperature. The maximum 
mass of the star is higher for higher entropies because the thermal effects make the 
EOS stiffer. 

Trapped neutrinos with a chemical potential /i^ are included by fixing the lepton 
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Fig. 2. Star mass as a function of radius for different configurations 



number defined as Yi = {pe + Pv^)/ Pb- In consequence there will be a large number 
of neutrinos in the star but also an increased electron density. Therefore, demanding 
charge neutrality, the proton density increases. For the configurations with 34 ^ 0, 
the high proton fraction delays the appearance of the proton chiral partner and 
for having less degrees of freedom the EOS becomes stiffer. For the configurations 
with 54 = 0, the chiral partners appears smoothly not having much effect, so the 
only effect in this case is that the high proton fraction makes the star more isospin 
symmetric and thus the Fermi energy smaller and the EOS softer. 

These two features can be put together to describe the evolution of the star. 
After the supernova explosion, the star is still warm so the entropy per baryon is 
fixed to S' = 2 (the temperature increases from at the edge up to 45 MeV in the 
center). The star still contains a high abundance of neutrinos that were trapped 
during the explosion so the lepton number is fixed to Yi = 0.4. After 10 to 20 
seconds the neutrinos can escape and beta equilibrium is established. After about 
one minute the temperature of the star has dropped below 1 MeV and the entropy 
per baryon can be considered zero. As it can be seen in Fig[3] and FigHl the effect 
of entropy and lepton number together is complicated. In the cases with finite (74 
(P2 and P4) the star's maximum mass decreases with time while in the case with 
54 equal to zero (PI and P3), the star's maximum mass increases with time. 

The problem is that this calculation does not take into account that the baryon 
number cannot change with time, otherwise the star would collapse into a black 
holJ^. If the baryon number is fixed the maximum masses for every configuration 
decrease with time. This maximum masses with fixed baryon numbers are repre- 
sented with dots in FigO and FigH) 

5. Rotation 

The maximum frequency with which a star can rotate without starting to expel 
matter on the equator, named Kepler frequency, has been determined by including 
monopole and quadrupole corrections to the metric due to the rotation and solving 
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Fig. 3. Star mass as a function of radius for different stages of the cooling for the configurations 
with Mjv_ = 1200 MeV 
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Fig. 4. Star mass as a function of radius for different stages of the cooling for the configurations 
with Mjv_ = 1500 MeV 

the self-consistency equation for as it was derived in . The rotation of the star 
generates a modification of the metric and the higher the rotational frequency, the 
higher is the mass and radius of the stai^l. For the different sets of parameters the 
increase in the maximum masses of the stars from v = to ~ vk fixing the baryon 
number to the value for zero frequency is of smaller than 5% (Fig. [5]), different from 
the 15% mass increase for the case that the baryon number is not fixed. But this 
situation can be identified as the spin down of a cold star with a certain baryon 
number that continues until it emits all its energy and stops rotating. At this point 
the star is considered "dead" . 



6. Conclusion 

With increasing density .i.e. toward the center of the star, chiral partners appear, 
reaching a point where they exist with about the same densities as the corresponding 
nucleons. The decrease in the scalar condensate signals the restoration of the chiral 
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Fig. 5. Star mass as a function of rotational frequency for different configurations 

phase. This transition is a cross over for any of the studied configurations due to 
the requirements of beta equihbrium and charge neutrahty. 

The maximum mass of the star is higher when the couphng constant 174 is set 
to zero (PI and P3) ahhough these values have to be constrained to a fixed baryon 
number during the different stages of the evolution - the constraint given by the 
warm case with trapped neutrinos for 174 = and for the cold beta equilibrium case 
for 34 7^ 0. In this way, the maximum mass and radius of the star decrease with 
time. A separate analysis shows that the mass and radius of the star increase when 
rotation is included. A study of cooling effects in the first seconds of the neutron 
star life together with rotation and angular momentum conservation is currently 
under investigation. 

References 

1. P. Papazoglou, D. Zschiesche, S. Schramm, J. Schaffner-Bielich, H. Stoecker and 
W. Greiner, Phys. Rev. C 59 (1999) 411.; S. Schramm, Phys. Lett. B 560 (2003) 
164; S. Schramm, Phys. Rev. C 66 (2002) 064310. 

2. C. DeTar and T. Kunihiro, Phys. Rev. D 39, 2805 (1989). 

3. D. Jido, Y . Nemoto, M. Oka and A. Hosaka, Nucl. Phys. A 671, 471 (2000) 
[arXiv:hep-ph/9805306 . 

4. D. Zschiesche, L. Tolos, J. Schaffner-Biehch and R. D. Pisarski, Phys. Rev. C 75 (2007) 
055202 arXiv:nucl-th /0608044_ . 

5. V. Dexheimer, S. Schramm and D. Zschiesche. larXiv:0710.4192l [nucl-th]. 

6. G. Baym, C. Pethick and P. Sutherland, Astrophys. J. 170, 299 (1971). 

7. T. Takatsuka, Nucl. Phys. A 588, 365 (1995). 

8. N. K. Glendenning and F. Weber, Phys. Rev. D 50, 3836 (1994)^ 

9. S. Schramm and D. Zschiesche, J. Phys. G 29, 531 (2003) arXiv:nucl-th/0204075] . 



